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OW  DISCRETE  VORTEX  SYSTEM  OF  RING  OF  FINITE  SPAN 

N.  F.  VoroMyev 

Here  we  investigate  the  problem  of  an  inviscid,  incoi pressible 
flow  past  the  lifting  surface  of  a wing.  The  wing  surface  itself  is 
replaced  by  vortex  surface  s,  while  the  sheet  of  vortices  flowing 
from  the  trailing  edge  and,  in  the  general  case,  fron  the  lateral  and 
leading  edges  of  the  wing,  is  represented  as  vortex  surface  £.  This 
surface  consists  of  vortices  whose  axes  in  the  case  of  steady  notion 
are  directed  along  the  flow  line.  Vortex  density  p*  on  surfaces  S and 
^are  determined  from  the  condition  of  nonpassage  (nonpenetration). 
Moreover,  the  conditions  of  shed  on  the  edges  of  the  wing  aust  be 
aet.  The  presence  of  a vortex  sheet  flowing  froa  the  edges  of  the 
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winq  makes  it  possible  to  assure  the  condition  of  velocity  limitation 
on  the  edges  from  which  the  vortex  sheet  flows  [1-3].  In  the  case  of 
a linear  3ystem,  when  the  vortex  sheet  flows  only  from  the  trailing 
edge,  on  the  leading  and  lateral  edges  of  the  wing,  as  we  know,  the 
speed  of  the  inflow  (penetration),  determined  within  the  framework  of 
an  ideal  fluid,  is  infinitely  great.  In  the  nonlinear  case,  where  the 
shape  of  the  surface  is  unknown,  solution  to  integral  equations  for 
wings  of  complex  plan  shape  is  difficult  [1,  2], 

There  exists  a method  for  calculating  the  aerodynamic 
characteristics  of  a flat  wing  of  arbitrary  plan  shape,  where  the 
vortex  layer  simulating  the  wing  surface  is  replaced  by  a system  of 
discrete  vortices,  the  intensity  of  which  is  determined  from  the 
conditions  of  nonpassage  (nonpenetration)  [4].  The  vortex  sheet 
outside  of  the  wing  is  also  simulated  by  discrete  vortex  lines,  which 
represent  a continuation  of  vortices  on  the  surface  of  the  wing 
itself.  Each  of  the  vortex  lines  outside  of  the  wing  consists  of 
rectilinear  segments  which  take  the  direction  of  velocity  at  the 
corresponding  point  in  space.  The  position  of  the  vortex  lines 
outside  of  the  wing  is  determined  by  the  method  of  successive 
approximations  in  the  calculation  process.  The  solution  to  the 
problem  of  the  flow  past  a wing  of  finite  span  according  to  the 
system  of  discrete  vortices  is  reduced  to  solving  a system  of 
algebraic  equations.  This  method  can  be  conveniently  used  on  the 
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computer  for  wings  of  arbitrary  plan  shape.  The  wing  can  also  be 
represented  by  a system  of  discrete  vortices  in  the  case  where  it  is 
cambered.  In  the  liscrete  systea  of  a wing  of  finite  span  the  problem 
of  conversion  in  the  case  of  an  increase  in  the  number  of  discrete 
vortices  replacing  t.he  wing  and  the  problem  of  satisfying  the  shed 
conditions  on  the  edges  of  the  wing  reaain  open. 

In  the  present,  study  we  show  that  with  proper  selection  of  the 
discrete  vortices  which  replace  the  wing  surface  and  points  at  which 
the  conditions  of  nonpassage  are  satisfied,  when  the  number  of 
vortices  is  increased,  the  algebraic  suns  used  to  represent  the 
velocity  induced  on  the  wing  surface  by  discrete  vortices  will  be 
transformed  into  integrals  whose  convergence  can  be  proven,  while  the 
introduction  of  additional  vortices  near  the  edqes  will  assure  that 
the  shed  condition  is  met. 

Generally  a wing  of  arbitrary  plan  shape  is  a certain  smooth 
surface  S,  which  must  be  covered  by  an  orthogonal  grid  of  curvilinear 
coordinates  related  to  the  wing  surface.  The  coordinate  system  is 
selected  such  that  the  line  £ = const  connects  the  leading  and 
trailing  edges  of  the  wing  (Fig.  1).  Applied  to  the  surface  is  the 
discrete  coordinate  grid  £ = coast,  C * const,  which  breaks  the  wing 
down  into  rectangles  with  sides  2Ay.  one  of  the  coordinate  lines 
passes  through  point  ^(6*.  Cj)  mad  is  the  coordinate  line  of  the  grid 
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5“li.  Til*  two  other  coordinate  lines  of  the  grid,  which  are  the 
closest  to  point  N(luls)  and  nr*  orthogonal  in  relation  to  lin*  i—lt, 
are  selected  such  that  point  ti)  ropr*s*nts  the  aiddl*  of  the  side 

of  the  coordinate  grid.  This  is  line  — At*.  Ci+A£j. 

Selection  of  the  direction  of  vortex  lines  which  replace  the 
wing  can  be  arbitrary.  These  vortices  are  not  subject  to  the  laws  of 
behavior  of  a vortex  line  in  the  flow  of  an  ideal  fluid.  It  is 
assuaed  that  the  segaent.  of  the  coordinate  line  whose  aiddle  is  point 
N(lu  5j)  is  th*  segnent  of  a vortex  line  of  constant  intensity.  Points 
Nfa,  tj),  through  which  the  n -shaped  vortices  pass,  can  be  nuabered  by 
row  and  colnnn:  every  point  has  a nnnber  (a,  n) . At  points 
(It,  l)  — Aw).  (£*,  A&)  the  vortex  li  ne  has  a break  and  continues  along 

coordinate  lines  tj— At*  t^Ci+ACj  to  the  trailing  edge  of  the  wing. 

Beyond  the  ^-shaped  vortex  of  intensity  AI\j  on  the  ning  on 
coordinate  line  segnent  t~=fcH-Afr— U+i,  the  niddle  of  which  is  point 
trf-Afc.  Cj-,  ve  find  the  following  vortex  line  of  constant  intensity 
AIWi,».  At  points  (S<+A£<,  t.(— Af/),  ($<+A|<,  t<+At<)  this  vortex  line  also  has 
a break  and  continues  along  coordinate  lines  A C*  C~trl-Atj  to  the 

trailing  edge  of  the  wing. 

Thus,  the  entire  wing  S is  covered  by  a systea  of  discrete 
rectangular  n-shaped  vortices,  which  are  related  to  the  wing.  The 
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n -shaped  vortices  arranged  on  the  wing  S,  and  also  a certain 
additional  number  of  discrete  vortices,  which  can  be  introduced  later 
for  meeting  the  conditions  of  shed  on  the  edges  of  the  wing,  descend 
froa  the  wing  and  continue  outside  of  the  wing  to  infinity, 
simulating  vortex  sheet  2.  The  vortex  lines  representing  vortex  sheet 
2 consist  of  finite  rectilinear  segments  which  lie  in  the  direction 
of  the  velocity  at  corresponding  points  outside  of  the  wing  [4].  The 
intensity  of  the  vortices  is  determined  from  the  conditions  of 
nonpassage  on  the  wing  surface  and  the  shed  conditions  on  the  edges 
of  the  wing.  Hhen  the  cell  dimensions  are  decreased,  in  the 
expression  for  the  velocity  induced  at  points  on  the  wing  surface,  we 
get  a peculiarity  related  to  sections  of  attached  vortices.  Belated 
to  the  sections  of  free  vortices  representing  vortex  sheet  2 is  the 
regular  part  of  the  expression  for  velocity  induced  at  points  on  the 
wing  surface  S.  for  proof  of  convergence  of  the  process  for  reduced 
cell  values  of  the  coordinate  grid  we  must  select  point  positions  on 
the  wing  surface  at  which  the  conditions  of  nonpassage  arm  satisfied. 
As  such  we  select  the  points  M (x,  x)  with  coordinates  *“=  •i+'y'i 
which  lie  in  the  geometrical  centers  of  the  coordinate  cells  (see 
Fig.  1). 

Henceforth  considerations  related  to  passage  to  the  limit  with  a 


decrease  in  the  coordinate  cell  A$(,  2ACj  without  limitation  of 
generality  will  be  done  for  a flat  wing  of  arbitrary  plan  shape 
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Everything  that  has  been  said  above  also  applies  to  arbitrary  smooth 
surfaces  without  curvature  separation  lines.  In  the  case  of  a flat 
wing  the  orthogonal  coordinate  systen  is  rectilinear  on  wings,  and, 
consequently,  each  n -shaped  vortex  related  to  the  wing  consist  of 
three  rectilinear  segment*. 


The  velocity  induced  by  elenentary  vortex  dl  of  intensity  T at  a 
certain  point  separated  by  distance  r fron  the  niddl®  of  the 
elenentary  vortex  is  determined  by  the  Biot-Savart  formula 

j(7 r Id/  y r] 

dV~irr-  r* 

Accord ing  to  this  formula,  for  points  on  a flat  wing  the  velocities 
induced  by  vortex  lines  lying  in  the  sane  wing  pl»ne  are  directed 
along  the  nornal  to  the  surface.  Here  th*  value  of  the  velocity 
induced  at  point  H(x,  z)  on  the  wing  by  a H -shaped  vortex  of 
intensity  ATg  passing  through  point  Wit.  Cj)  and  consisting  of 
rectilinear  segseats  of  finite  length  2AC*  5<  — *•(£;  — AC*).  h— x*(Cj+ACj)# 
where  6— *•  (0  is  the  equation  for  the  trailing  edge  of  the  wing,  can 
be  represented  in  the  fora  of 


AVi ) — — r,  5c  C>  -r  AC,)  — F ( x . 5c  C> — AC,)) 

u in®  (St-  C i *)• 


where 
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Fix,  z,  l,  ;>  - 


1 I ) U-E)»+(«-i 


l*-»,  (01 

VI*-*,  OP +<«-.') 


(*.  *.  . * + e>)  — F (x,  Z,  Si . * — *5>)1 , 

tJ*»/  (It.  &>  = *). 


where 


Ari=Arl;  (£,.  z),  Afc-v,. 


Point  H(x,  z)  is  selected  in  the  center  of  the  coordinate  cell, 
so  that  ia  the  case  of  a finite  nuaber  of  discrete  vortices  Sis*1*, 
while  5j —z  only  for  a certain  column  of  coordinate  cells.  The 
velocity  induced  it  point  n (x,  z)  by  all  discrete  n -shaped  vortices 
associated  with  the  wing  S is  represented  by  the  sum  of  velocities 
induced  by  each  of  the  n -shaped  vottices.  If  we  perform  summation  in 
a fixed  row,  in  this  case  discarding  the  term  which  correspond  to 
value  £ = z,  and  then  suaaation  for  all  rows,  the  velocity  at  point 
(x,  z)  can  be  represented  in  the  fora  of 


v = ~ lir (2 2 p M h-'i-Kf)  x 

1 » B 

y 2A£/A|t  + 2p(gi,  z)(F(x,z,  * + ej) - F(x,  z,  git  z-e,)l  A6,|, 

where  the  intensity  ATy  of  the  n-shaped  vortex  passing  through  point 
5*.  li  (U.  M is  represented  in  the  form  of  Arij«p(E/,  {j)  Air,  here  i*  the 
distance  along  axis  f from  point  £<•  U)  of  this  vortex  to  point 

&(W+>>  W#  th  rough  which  the  following  fl -shaped  vortex  passes. 
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In  foraula  (3)  sunaation  by  coluan  is  excluded  in  the  doable  sub, 
where  C"“*.  Sunaation  of  this  coluan  is  done  separately. 


For  an  infinite  increase  in  the  nunber  of  discrete  vortices, 
•hen  first  2A£|-*0,  and  then  ACj-eO,  the  velocity  induced  at  point  (x, 
by  all  attached  vortices  can  be  represented  in  the  fora  of 


V = - 


"iTT  ( XT  P(E.  Qft(x,  *,tOdCdi+  J p(5,  Z)IF(X,  *,  1,  *-fe)- 

l*— *•  *u  (*) 


— F(x.  z,  t,  i — e))ig 


(4) 


where  S-2£  is  the  area  of  the  wing  with  the  excluded  flat  width  of  2£ 
near  point  »**•  £ =*»(C).  ? = **(C)  represents  the  equation  of  the  trailing 
and  leading  edges  of  the  wing,  respectively,  while  the  derivative  of 
the  function  P(x,  z,  5,  £)  takes  the  form  of 


Fi(x,z,  5,  0 = 


L 


<*  ~ [V(T-  £)'+<*-{)' 

_ i« -«,(£»»+»  i*  - x,  (oi*  (*  - - *;  o (« _ :)i 

<i-*-xJ(j)p+  (i~~£)7ys! 


Function  Fc(*»*  -1.0  exists  in  the  range  of  S-2  £.  In  the  range  of 

S-2t  within  the  internal  integral  of  the  double  integral  of  foraula 
(«)  we  caa  integrate  by  parts.  After  setting  z > 0,  without  Uniting 
the  generality  we  get 

I f rx  V°> 

^ = l pI5,*,(01x 

*s  <•*  *n  (x— •) 

X Fix. z,t 2, ($)]<£  + I p(5,0)F  (x,  2,1,0)  ft  + f p [|.  znp  (l)]  : 
o o 

*j(0> 

X F u,  2,  t 2„P  (1)1  *5  + I pit  2up  (5)1  F lx,  Z,  t Znp  (5)]  di  - 

x (x- •) 
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*a  U— •)  JT,  (0) 

- J P(iO)f(x,*,i^)^-  f p(l0)f(x.*,5.0)d$  + 

* JT.  (*-*) 


«,(*—•)  «,(*+*>  >C| 

--  ]■  p$.z.-e)F  (x,2,l,z  — C)dt—  \ p (1,0)  F(x,  2,1.0)  dt-t  1 1 

*u  <*•—*>  *n  («-  *) 

*,<»+«)  x,(jf«) 

+ J p IE.  *op  (6)1  F fx,  *,S^np(5)1  — f p(i,z-!- e)F(x,2,l,z -ri)dl^ 

*m  w •>  *■  (*+•) 


**  <»>  1 

+ f p(i*)l/f(x,*,ii  + e)  — F(x,i,t*  — e)ld5  • 

*«<*)  ' 

For  proof  of  ths  existence  of  velocity  on  the  wing  determined  by 

foreala  (S),  we  eust  Bake  assumptions  on  the  shape  of  the  wing 

contour  and  on  the  nature  of  vortex  density  oe  the  wieg.  For  contour 

L we  assune  the  continuity  of  equations  £“*»(£)  of  the  leading 

and  trailing  edges  of  the  contour  in  the  range  of  a < C < b (can  be 

the  lateral  edges,  parallel  to  axis  C when  £ = a,  £ = b)  . He  also 

assuae  that  on  tha  wing,  including  the  edges  of  the  wing,  the  value 

of  vortex  density  p(£,  £)  and  derivative  P<  (£,£),  from  whose  value  we 

determine  the  naxiaal  density  of  discrete  vortex  lines  on  wing  S 

coinciding  with  the  direction  of  axes  £,  satisfy  the  Holder  boundary 

condition. 


Wow  let  us  prove  the  existence  of  velocity  on  a wing  determined 
by  foraula  (5)  for  internal  wing  points.  The  double  integral  in 
foraula  (5)  can  be  represented  in  the  fora  of 


where  expression 


£)•(*.  *.LCMU5 


V (T-  +Ji -T? V «, (W -E (« ~ t)‘  - (x - S) [x - x, (0 1 


(6) 
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represents  the  continuous  function  for  winq  points  (x,  z)  which  do 
not  lie  on  the  trailing  edge  of  the  wing.  This  integral  represents 
the  main  value  of  the  Cauchy  type  iterated  integral,  which  exists  for 
points  which  do  not  lie  on  the  wing  contour  and  which  because  of  this 
represents  the  internal  points  for  each  of  the  iterated  integrals 

m. 


After  adding  the  single  integrals  in  which  we  have  the 
expression  p(F,  0)  F(x,  z,  F,  0)  dF  under  the  integral  sign,  there 
remain  two  teras,  which  on  the  basis  of  the  main-value  theorem  can  be 
represented  in  the  fora  of 


*n(*  *••)  at,  (/— «) 

J P&0)F(x,z,lO)dt  + J p(5,0)F(x,z,5.0)rf5  = 

xn  <*— *)  i,  (»+*) 

= p (*).  0)  F \x,  z,  x„  (z),  0)  2e  + P lx,  (z),  0)  F lx,  z,  x3  (z),  01 2e, 

F \x , z,'xB't  (»),  01  = 1 j Of  '"•»  1 

* l * *■.»  W |/|*  — (0) J*  -H »•  | 


is  the  function  United  to  points  which  do  not  lie  on  the  edges  of 
the  wing  (value  z > 0)  . Under  the  above  assumptions  on  the  finiteness 
of  the  values  of  vortex  density  on  the  wing  edges,  the  value  of  each 
of  these  terns  when  e-*-0  rswsrts  to  zero. 
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The  sub  of  single  integrals  which  under  the  integral  sign 
contain  function 


f <«.». i.± .)  ' (VH3H2  ) , 

' (T  •)}  *-l  V\x~  *,(*  ± •)!■  + «*  I 


where 


_ ± «)  . 
(i*5 

lim//(x,2,&z  + e)  =0. 


also  reverts  to  zero  when  Let  us  show  this  using  function  P(x, 

z,  I*  z ♦ £ ) as  an  exaaple. 

*,(1)  «,<*+«) 

f p(i*)f  (x,z,$,z  + e)dg  - f p(S,z  + t)F(x,z,i2.f  «)4» 

*«  (*)  *■  <»+*> 

*,  <•+«)  *B  <*+«) 

= — J [p(E,2  + e)  — p(S,2)J/r(jC,2,6,2  + t)dl+  J p(6,2)X 

rn  »+«>  *n  <*> 


«,  (*) 


V'+«>  r . 


XF(*,:,5,z  + e)4+  J p(iz)F(x,z,6,2  + e)d{  = f [epc(t*)  + 

«,  U-H)  V&+** 


*B  <*>  + «BW 


+ 0(8^1 (% i_  f p(5,2)//(x,2,6,2  + e)rf| 

*o<*> 

I 

j p(5,2)//(X,*,i2  + t)d5. 

*j  (l*  f 

Tn  the  first  tern  the  integrand,  because  of  the  property  of  function 
H,  reverts  to  zero  when  e-»-0#  while  the  last  two  terns  are 
representad  on  the  basis  of  the  aeaa-value  theoren  in  the  fora  of 

— P |*o  (?)  -r  -J-  x a (*)•  *|  H Jx,  2,  Xn  (z)  + -j-  XB  (z),  2 -4  ej  X,  (2)  — 

- P Jx,  (2)  -f  -J-  xi  (2),  zj  //  [x,  2,  X,  (2)  + -j-xi  (2),  2 + t j X,  (z) 


r 


one  - 1790 


PARE  12 


and,  on  the  strength  of  function  H,  also  revert  to  zero  when  e-*0. 

The  reaaining  single  integrals  in  the  right  part  of  fornula  (5) 
represent,  when  a-eO  the  contour  integral 

where  5 = f(f)  is  the  equation  of  the  contour  L,  integration  with 
respect  to  which  is  done  counterclockwise.  The  contour  integral  can 
he  written  in  the  forn  of 

f MU  (1)1  «• 

where  function  <D  is  detersined  by  dependence  (6).  For  points  (x,  z) 
which  do  not  lie  on  contour  L with  the  condition  of  continuity  of 
equations  for  the  leadinq  and  trailing  edges,  two  situations  can  be 
encountered  in  the  contour  integral:  1)  x = €,  z t f (f) , 2)  x / C,  z~ 
= f (§•) . The  case  where  x = 5,  z = f(€)  for  internal  points  (x,  z)  of 
a wing  cannot  exist. 

A F,  (6)  d\ 

In  the  first  case  the  integral  jr  /—  \~ , shore 
F,  ({)  = P It represents  a continuous  function,  where  Pt  (x)  = 
pf  x,  f (x)  ] is  a Cauchy  type  integral  and  exist  in  the  sense  of  the 
•ain  value. 


In  the  second  case  the  integral  tm  shore 
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fi(0  = la  a continaoas  function.  Here  P,(C°)  = 0 for 

n -shaped  vortices  which  and  on  tha  trailing  edge  is  reduced  to  the 
fora  of  the  Cauchy  type  integral  by  substitution  of  variables  f (5)  = 


where  when  f ->  t — ^ z.  This  integral  exists  in  the  sense  of  the 
main  value  when  f*  (5°)  f 0.  Por  contour  L with  the  condition  of 
continuity  of  the  equations  of  the  leading  and  trailing  edges  value 
f* (?)  = 0 can  occur  only  for  end  points  on  the  contour  £ * a,  £ = b, 
which  cannot  be  points  ? = c®  for  points  (x,  z) , which  do  not  belong 
to  contour  I..  If  there  are  lateral  edges  parallel  to  axis  ? when  £ = 
a,  £ * b,  where  f*  (?)  = 0,  then  neither  can  the  points  on  the  lateral 
edqes  be  points  ? 3 ?°  for  points  (x,  z)  which  do  not  belong  to  the 
wing  contour.  This  Beans  that  in  the  second  case  the  contour  integral 
exists  in  the  sense  of  the  sain  value. 

Thus,  the  velocity  induced  by  the  n -shaped  vortices  on  wing  S 
with  transition  to  the  lisit  froa  the  discrete  system  to  the  systea 
of  the  vortex  surface,  is  deterained  for  internal  wing  points  by  the 
foraula 


v = lH*r  ^&QFl*'*'t'VdW+  (7) 


i 
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The  integrals  which  deternine  velocity  are  integrals  of  the 
Cauchy  type,  and  for  internal  wing  points  they  exist  in  the  sense  of 
the  aain  value.  Essential  for  representing  velocity  in  the  sense  of 
the  aain  value  of  Cauchy  type  integrals  (including  iterated 
integrals)  is  the  breakdown  of  the  integration  range  into  parts 
within  the  neighborhood  of  a certain  point:  0 ^ 5 ^ i - 6,  x ♦ 5V<  ( 
^ x3(0)  and  a^C^z  - £,  z ♦ £ ( v<  b.  The  selection  of  the 

point  fl(x,  z)  in  the  center  of  the  coordinate  systea,  carried  out  in 
this  work  for  the  discrete  case,  assures  convergence  of  the  integrals 
in  the  sense  of  the  main  value  with  transition  to  the  liait. 

The  end  (nonzero)  density  values  of  the  vortex  lines  on  portions 
of  the  wing  contour  where  the  vortex  sheet  does  not  descend  (flow 
off),  give  us  infinite  velocity  values  within  the  franework  of  an 
ideal  fluid.  On  wing  contour  L the  conditions  for  the  existence  of 
the  integrals  through  which  the  velocity  induced  by  the  vortices  is 
expressed  will  be  net  in  the  case  where  contour  L is  a line  which 
lies  entirely  within  the  vortex  surface  S ♦ 6S.  This  neans  that  the 
vortex  surface  of  the  wing  aust  continue  unbroken  beyond  wing  S.  Then 
contour  L becones  a line  whose  points  are  the  internal  points  of  the 
vortex  surface  S ♦ 6S  which  lie  within  the  contour  L ♦ 6L.  For  the 
internal  points  of  the  surface  the  existence  of  velocity  in  the  sense 
of  the  aain  value  is  proven.  Here  the  velocities  at  points  on  contour 
L are  deternined  by  foraula  (7),  where  integration  is  done  with 
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respect  to  the  surface  S ♦ 6S  and  the  contour  L ♦ 6L.  The  vortex 
sheet  lyinq  within  the  contour  L ♦ 6L,  whose  shape  in  the  nonlinear 
case  is  known  in  advance,  will,  as  already  nentioned,  give  us  the 
regular  velocity  component  for  points  on  the  wing  S and  its  edges  L. 

In  the  studied  case,  where  the  wing  is  a lifting  surface  without 
thickness,  for  the  existence  of  finite  velocity  on  the  edges  we  nust 
iepose  the  condition  of  the  seooth  joining  of  the  wing  surface  s and 
the  vortex  sheet  2 which  flows  off  of  it  and  the  condition  of 
continuous  transition  of  the  vortex  density  of  these  surfaces  on 
their  boundary  - contour  L.  On  contour  L satisfaction  of  the 
condition  of  finite  velocity  does  not  generally  require  that  the 
vortex  density  on  the  edge  of  the  wing  revert  to  zero,  provided  the 
vortex  sheet  flows  frow  it.  The  vanishing  of  the  density  of  the 
vortex  lines  which  coincide  with  the  direction  of  axis  5 on  the 
trailing  edge  of  the  wing  in  the  case  of  the  linear  systen  is  caused 
by  the  fora  of  the  vortex  sheet  beyond  the  wing  and  coaes  froa  the 
condition  of  shed  froa  the  edges  of  the  wing  foraulated  above.  In  the 
linear  systea,  where  the  vortex  lines  beyond  the  wing  take  the 
direction  of  velocity  at  infinity,  on  surface  2 there  are  no  vortex 
line  coaponents  which  are  perpendicular  to  axis  f (axis  6 on  the  wing 
as  the  direction  of  velocity  at  infinity)  , and  froa  the  condition  of 
the  continuous  transition  of  the  vortex  surface  S into  surface  I it 
follows  that  on  the  trailing  edge  of  the  wing  S the  intensity  of  the 
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vortex  density  of  the  last  of  the  H -shaped  vortices  should  be  equal 
to  zero.  Por  the  vortex  line  (wing  of  infinite  span)  fron  the 
condition  of  finite  velocity  we  also  get  the  values  for  zero  vortex 
density  on  the  trailing  edge  of  the  vortex  segnent.  In  the  case  of 
the  end  point  of  a vortex  line  the  condition  of  finite  velocity, 
detemined  by  the  Cauchy  type  integral,  can  be  satisfied  only  when 
vortex  density  vanishes  at  the  end  of  the  line  [5].  However,  in  the 
case  of  a vortex  shedding  fron  the  edges  of  a wing  of  finite  span, 
the  edges  of  the  wing  do  not  represent  the  end  points  of  the  vortex 
lines,  and  vortex  density  on  the  shed  line  is  generally  not  equal  to 
zero. 

In  the  discrete  systea  we  can  iapose  a distribution  of  vortex 
lines  which  at  the  transition  to  the  liait  froa  the  discrete  systea 
to  the  systea  of  the  vortex  surface  would  cause  a break  in  the 
transition  of  vortex  density  on  the  wing  contour  L.  This  is  assured 
if  a point  which  in  the  discrete  systea  is  considered  to  be  a wing 
edge  point  is  liaited  on  the  vortex  sheet  side  by  a vortex  line  of 
the  saae  intensity  as  on  the  inner  side  of  the  wing  (Pig.  2).  Since 
on  wing  S the  vortex  lines  parallel  to  axis  £ and  axis  ( have  two 
directions,  for  the  selected  systea  of  rectangular  n -shaped 
vortices,  the  liaiting  density  value  of  discrete  vorticea  on  the  wing 
In  the  direction  of  axis  C equals  p(£  = const,  C)  * while  the  liaiting 
density  value  of  liscrete  vortices  on  the  wing  in  the  direction  of 

] 
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axis  C equals  const)  dp , for  each  point  on  the  edge  at 

which  the  condition  of  nonpassage  is  satisfied  there  are  generally 
two  additional  vortices.  On  the  side  of  the  vortex  sheet  one  vortex 
runs  parallel  to  axis  p,  the  other  - parallel  to  axis  C. 


The  wing  is  broken  down  by  a discrete  grid  of  coordinate  lines 
into  a coordinate  cells.  In  the  center  of  each  of  these  at  point  (a, 
n)  the  condition  of  nonpassage  is  satisfied.  Here,  on  contour 
segaents  which  do  not  coincide  with  the  direction  of  the  coordinate 
lines,  the  contour  is  replaced  by  a broken  line,  which  is  located 
outside  of  th*  wing  (the  wing  surface  is  taken  with  an  excess)  . The 
rectangular  vortex  lines  which  are  bound  to  the  wing  and 

whose  intensity  la  deterained  froa  the  conditions  of  nonpassage,  are 
shown  in  the  left  side  of  Pig.  2 as  continuous  lines.  The  arrows 
indicate  the  positive  direction  selected.  The  figure  also  shows  the 
part  of  these  lines  which  lies  in  a single  cell.  Here  each  of  the 
bound  vortex  lines  begins  and  ends  on  the  trailing  edgp  of  the  wing. 
Continuation  of  these  vortex  lines  outside  of  the  wing,  beginning 
with  the  leading  edge,  coincides  with  the  direction  of  the  flow 
velocity.  On  the  left  side  of  the  figure  their  continuation,  now  as 
free  (vortex  lines),  is  shown  by  solid  lines  which  begin  at  the 
corresponding  points  on  the  trailing  edge.  The  boundary  points  at 
which  the  conditions  of  nonpassage  are  satisfied  and  which  for  the 
discrete  systea  represent  the  contour  points,  are  aarked  by  x's.  At 
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these  points  the  shed  conditions  nust  be  aet,  i.e.#  conditions 
iaposed  which  assure  continuous  transition  of  the  vortex  surface  of 
the  wing  S into  the  vortex  surface  of  the  sheet  2. 

For  a trailing  edge  which  is  not  parallel  to  coordinate  axis  C# 
point  (a,  1)  which  lies  in  the  center  of  rectangle  is 

Uniting.  Here  the  boundary  of  the  wing  S and  of  the  vortex  sheet  2 
are  the  sides  and  to  assure  an  unbroken  vortex 

surface  in  the  direction  of  axis  ? (vortex  lines  coincide  with 
direction  of  axis  C)  along  the  side  Dm.,Cm,\  , vortex  line  of  intensity 
AIVi  is  introduced,  where  represents  the  intensity  of  the 

H-shaped  vortex  line  bound  to  the  cell  whose  center  is  the  point  (a, 
1)  . At  points  Dm.\,  Cm, i the  vortex  line  breaks  and  thereafter  behaves 
as  a free  vortex  line.  To  assure  the  continuity  of  the  vortex  surface 
in  the  direction  of  axis  C along  side  t ««  introduce  the  vortex 

line  of  intensity  2AlVi  — At  points  the  introduced  vortex 

line  breaks  and  thereafter  behaves  as  a free  vortex  line.  In  Fig.  2 
the  introduced  vortices  are  shown  as  dashed  lines.  There  selected 
positive  direction  is  indicated  by  arrows. 

For  a trailing  edge  parallel  to  axis  C point  (n,  n)  , which  lies 
in  the  center  of  the  rectangle  Bm,*Fm,HD  is  the  boundary  point. 

Here  the  boundary  of  the  wing  S and  of  the  vortex  sheet  2 is  the  side 
Bm,nCm,nm  fo  assure  the  continuity  of  the  vortex  surface  along  side 
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, we  introduce  a vortex  line  of  intensity  AF*.*.  it  points 
the  introduced  line  breaks,  and  thereafter  behaves  as  a free 

vortex  line. 

[ 

For  a lateral  edge  parallel  to  axis  c.  the  point  (a,  1),  which 
lies  in  the  center  of  rectangle  Bm_iCm.tDm,x  Fm,i , is  the  boundary  point. 
Here  the  boundary  of  vinq  S and  of  the  vortex  sheet  ~ is  side  - 
Along  side  we  introduce  a vortex  line  of  intensity  ‘2Arm,t — AF,,.*, 

which  continues  to  the  end  of  the  lateral  edge,  whero  it  sheds  as  a 
free  vortex  line  together  with  all  attached  vortex  lines  coinciding 
with  the  lateral  edge.  At  point  Bm,i  the  vortex  line  behaves  as  a free 
vortex  line.  The  selected  positive  direction  of  the  bound  and 
introduced  vortex  lines  is  indicated  by  arrows. 

For  a leading  edge  which  is  not  parallel  to  axis  C the  point  (n, 
0),  which  lies  in  the  center  of  rectangle  D*,o  #*,u  Fm.u  C*,* , is  the 
boundary  point.  Here  the  boundary  of  the  wing  S and  of  the  vortex 
sheet  - consists  of  the  sides  Dm.oB and  0m.oC„,,o„  Along  side 
we  introduce  a vortex  line  of  intensity 
AFm+i,i  — Arm,t (AFm . ,.t  — AFm,i  -I-  AF,n>0)y  which  fros  point  Dm.t  behaves  as  a free 
vortex  line,  and  froa  point  Bmt(Cm, o>  continues  along  the  wing  parallel 
to  axis  ? to  the  trailing  edge  of  the  slag,  where  it  sheds  as  a free 
vortex  line  together  with  the  corresponding  attached  vortices.  The 
positive  direction  of  the  introduced  and  attached  vortex  line  is 
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indicated  by  arrows.  The  different  signs  for  the  value  of  intensity 
in  the  case  of  corresponding  arrows,  and  also  identical  signs  in  the 
case  of  arrows  of  opposite  directions,  indicate  the  mutual 
obliteration  of  the  effect  of  the  vortex  lines.  As  we  see  in  Fig.  2, 
the  intensities  of  the  introduced  vortex  lines  which  assure  the 
continuous  transition  of  vortex  surface  S into  surface  2 in  the 
direction  of  axis  £ and  axis  5 for  point  (m,  0),  are  such  that  they 
Mutually  obliterate  the  effect  of  the  attached  vortex  bound  to  point 
<■#  0)  and  the  effect  of  the  sums  of  introduced  vortices  of  intensity 
Ar«.o.  in  the  discrete  wing  system  it  is  possible  to  not  immediately 
introduce  the  attached  vortex  of  intensity  Ar*,»  which  corresponds  to 
(a,  0),  but  to  select  the  intensities  of  the  introduced  vortices 
equal  to  AIV-u  and  — Ar-m,i  f respectively. 

Thus,  the  intensity  of  all  newly  introduced  vortices  on  the  wing 
edqes  which  assure  finite  velocity  at  the  points  on  the  edges,  is 
expressed  as  the  intensity  of  the  attached  H-shaped  vortices  plotted 
earlier,  whose  intensity  is  determined  from  the  conditions  of 
nonpassage  on  the  wing  surfaces. 
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